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SUMMARY 


\ 


A  survey  of  available  wind  shear  data  was  conducted  and  the  results  were  re¬ 
lated  to  current  design  practices  for  vertically  rising  missiles. 


the  responses  of  several  different  kinds  of  vehicles  were  obtained  by  means  of 
a  6-dimensional  trajectory  computer  program. 

The  maximum  wind  velocity  and  the  integrated  area  under  the  wind  profile 
(up  to  the  critical  altitude  at  which  the  maximum  wind  velocity  occurs)  were 
found  to  be  critical  parameters.  A  design  procedure  was  developed  for  the 
construction  of  a  design  response  diagram.  This  relates  the  vehicle  response 
to  the  maximum  wind  velocity,  its  probability  distribution,  and  the  probability 
distribution  of  the  integrated  area.  Since  no  statistics  in  the  integrated  area 
are  available,  a  method  of  estimation  was  also  developed,  involving  the  use 
of  four  design  wind  profiles  ^Ap^ewdtsr^C^^These  design  diagrams  are  useful 
during  all  phases  of  design,  as  well  as  pre-launch  checkout. 

A  detailed  discussion  is  presented  of  the  probability  concepts  involved  showing 
that,  in  most  cases,  the  probability  of  occurrence  caijxSnly  be  defined  by  "much 
less  than"  a  certain  percentage. 

Several  recommendations  for  futur^-sfudy  are  made.  These  include  reduction 
of  wind  sounding  data  to  obtertlfmte  grated  area  statistics  and  further  refine- 
menigjRl'ffie  design  procedure  for  winds  aloft. 

The  proposed  procedure  has  a  high  degree  of  flexibility  which  permits  optimi¬ 
zation  of  the  design  of  a  missile  system.  As  a  result  the  launch  probability 
can  be  optimized  with  respect  to  the  mission  of  the  vehicle  and  the  structural 
weight  traded  off  against  payload  capability  or  flight  performance. 


111 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


LMSC- A003108A 


CONTENTS 

Page 

SUMMARY  ili 

TABLES  vii 

ILLUSTRATIONS  vi 

SECTION  1  INTRODUCTION  1 

SECTION  2  WIND  VELOCITY  DATA  4 

SECTION  3  VEHICLE  RESPONSE  7 

Parameter  Study  8 

Design  Profiles  10 

Effect  of  Different  Vehicles  12 

Summary  of  Parametric  Relationships  12 

SECTION  4  DESIGN  PROCEDURE  14 

Development  14 

Detailed  Procedure  15 

Effect  of  Wind  Direction  16 

Pre-Launch  Checkout  Procedure  17 

SECTION  5  PROBABILITY  CONSIDERATIONS  19 

Wind  Profiles  19 

Design  Diagram  20 

SECTION  6  CONCLUSIONS  AND  RECOMMENDATIONS  23 

SECTION  7  REFERENCES  26 

APPENDIX  A  DEVELOPMENT  OF  DESIGN  WIND  PROFILE  57 

SERIES 


v 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


LMSC- A003108A 


ILLUSTRATIONS 


Figure 

Page 

1 

Annual  Winds  at  Patrick  AFB 

32 

2 

Summary  of  AN/GMD-2  Wind  Data 

33 

3 

Sissenwine  Wind  Profiles  for  1  Percent 

Probability  During  Winter 

34 

4 

Typical  Design  Wind  Profiles  for  1  Percent 
Probability  During  Winter 

35 

5 

Avidyne  Wind  Profiles  for  1  Percent 

Probability  During  Winter 

36 

6 

Marshall  SFC  Wind  Profiles  for  Saturn  Vehicle 

37 

7 

Tail  Wind  Response  for  Various  Wind  Profiles 
(Vehicle  A) 

38 

8 

Response  for  Wind  Profiles  from  Literature 
(Vehicle  A) 

39 

9 

Synthetic  Wind  Profiles  for  Parameter  Study 

40 

10 

Response  for  Parameter  Study  Profiles 
(Vehicle  A) 

41 

11 

Regression  Curve  of  Rigid  Body  Side  Wind 

Response  (Vehicle  A) 

42 

12 

Regression  Curve  of  Rigid  Body  Head  Wind 

Response  (Vehicle  A) 

43 

13 

Regression  Curve  of  Rigid  Body  Response  for 
Maximum  Wind  Velocity  of  300  fps  (Vehicle  A) 

44 

14 

Rigid  Body  Response  for  Design  Wind  Profiles 
(Vehicle  A) 

45 

15 

Rigid  Body  Response  for  Design  Wind  Profiles 
(Vehicle  B) 

46 

16 

Rigid  Body  Response  for  Design  Wind  Profiles 
(Vehicle  C) 

47 

17 

Design  Wind  Response  ^Vehicle  A) 

48 

18 

Design  Wind  Response  (Vehicle  B) 

49 

19 

Design  Wind  Response  (Vehicle  C) 

Vi 

50 

LOCKHEED  MISSILES  ft  SPACE  COMPANY 


LMSC-A003108A 


ILLUSTRATIONS  (Continued) 

Page 

20  Design  Diagram  for  Wind  Shear  Response  51 

(Vehicle  A) 

21  Design  Diagram  for  Wind  Shear  Response  52 

(Vehicle  B) 

22  Design  Diagram  for  Wind  Shear  Response  53 

(Vehicle  C) 

23  Shear  versus  Shear  Length  for  Wind  Profiles  with  54 

Maximum  Wind  Velocity  of  300  fps 

24  Shear  versus  Shear  Length  for  Wind  Profiles  with  55 

Maximum  Wind  Velocity  of  225  to  260  fps 

25  Parameters  for  Design  Wind  Profiles  D1  64 

26  Typical  Examples  of  Design  Wind  Profiles  D3  and  D4  65 

27  Parameters  for  Design  Wind  Profiles  D3  and  D4  66 

TABLES 

1  Summary  of  AN/GMD-2  Wind  Data  28 

2  Wind  Shear  Summary  29 

3  Maximum  Response  for  Various  Wind  Profiles  30 

(Vehicle  A) 

4  Synthetic  Wind  Profiles  for  Parameter  Study  31 

5  Parameters  for  Design  Profiles  D1  61 

6  Parameters  for  Design  Profiles  D2  62 

7  Parameters  for  Design  Profiles  D3  and  D4  63 


vii 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


LMSC-A003108A 


SECTION  I 
INTRODUCTION 


It  is  widely  recognized  that  wind  shear  and  gust  are  important  and  often  de¬ 
cisive  criteria  in  the  design  of  missile  systems.  This  unanimity  is  lacking, 
however,  when  an  answer  is  sought  to  tne  question:  "What  constitutes  adequate 
and  sufficient  criteria  for  wind  shear  and  gust  in  missile  design.  "  General 
agreement  is  that  such  criteria  must  have  a  probabilistic  basis  but  the 
interpretation  of  incomplete  statistical  data  is  the  cause  for  much  diversity  of 
opinion.  Principal  sources  of  this  confusion  are  listed  below: 

a.  The  extreme  variability  of  atmospheric  winds 

b.  The  lack  of  statistically  adequate  data,  not  only  for  the  American 
continent  but,  to  a  lesser  extent,  also  for  specific  locations 

c.  The  poor  reliability  of  much  of  the  data  collected  to  date,  princi¬ 
pally  due  to  instrument  errors  (Refs.  1  and  Z) 

d.  The  various  means  employed  to  reduce  the  raw  data  obtained  by 
soundings  (Refs.  Z  and  3) 

The  generally  accepted  procedure  for  evaluating  wind -Induced  structural  loads 
is  to  approximate  the  extreme  vehicle  response  using  wind  profiles  of  given 
estimated  probabilities.  A  wind  profile  is  a  graph  of  wind  velocity  as  a 
function  of  altitude,  where  wind  velocity  is  defined  by  the  method  used  for 
the  analysis  of  wind  sounding  data.  The  earliest  procedure  was  to  analyze 
the  radiosonde  data  statistically  for  horizontal  wind  velocity  only.  The  re¬ 
sulting  profiles  depict  average  wind  velocities  and  average -plus -n-number- 
standard -deviations  extreme  winds,  (Fig.  1).  This  process  smoothes  out  the 
wind  shears  and,  since  shears  are  critical  for  vertically  rising  vehicles,  it 
destroys  the  usefulness  of  the  data  for  missile  design.  The  more  recent 
approach  is  to  compute  the  shears  for  each  sounding  and  then  evaluate  the 
shears  statistically  (Refs.  1  and  2).  Since  only  shears  are  obtained,  the  pro¬ 
files  need  to  be  rounded  out  with  velocity  statistics  obtained  by  the  former 
procedure  (Refs.  7  and  8). 
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Several  factors  have  led  to  the  general  acceptance  of  wind  profiles  for  engi¬ 
neering  application: 

a.  Maximum  wind  velocity  and  shear  generally  occur  in  the  altitude 
interval  of  30,  000  to  45,  000  feet,  the  same  region  where  critical 
flight  conditions  usually  occur. 

b.  Wind  direction  remains  fairly  constant  with  altitude  when  the  wind 
velocity  is  well  above  average  (Ref.  8). 

c.  Strong  wind  velocities  and  strong  vertical  wind  shears  are 
associated  with  each  other,  at  least  in  the  levels  of  maximum 
wind  velocity  in  the  upper  troposphere  (Refs.  1  and  8). 

d.  Determination  of  the  vehicle  response  to  a  wind  profile  input  is 
rather  straightforward  and  the  effects  of  gust  and  elastic  body 
modes  may  be  obtained  by  superposition. 

e.  Sufficient  radiosonde  data  are  available  for  determining  wind 
velocities  and  shears  at  various  confidence  levels  with  a 
reasonable  degree  of  reliability. 

With  the  selection  of  the  procedure  utilizing  wind  profiles,  it  remains  to 
determine  specific  design  profiles  for  different  levels  of  probability.  Inspec¬ 
tion  of  the  wide  variety  of  one  percent  wind  profiles  published  in  the  literature 
(Refs.  1,  7,  8,  9,  10,  11,  12,  13,  14)  shows  that  there  is  general  agreement 
that: 

a.  Maximum  wind  velocity  and  shear  should  coincide; 

b.  Maximum  wind  velocity  generally  occurs  in  the  altitude  interval 
between  30,  000  to  45,  000  feet. 

c.  The  vehicle  should  be  analyzed  for  wind  shear  at  least  at  the  most 
critical  altitude  within  the  above  interval. 


In  addition,  reasonably  close  agreement  exists  on  the  magnitudes  of  the 
maximum  wind  velocity  and  shear.  On  the  other  hand,  areas  of  disagreement 
are: 

a.  Should  a  gust  be  superimposed  upon  the  shear  and  if  so,  what  should 
be  the  gust  velocity  and  gust  length. 

b.  How  does  the  wind  velocity  vary  for  altitudes  below  and  above  the 
shears. 
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c.  Should  the  dealgn  wind  profile  incorporate  aero  ground  wind  or 
statistical  ground  wind. 

Since  the  objective  of  this  study  is  the  development  of  wind  criteria  for  verti¬ 
cally  rising  vehicles  and  not  the  general  investigation  of  atmospheric  winds, 
the  determination  of  a  design  procedure  was  made  dependent  upon  the  response 
of  several  current  missiles  to  a  variety  of  wind  profiles.  A  study  of  critical 
response  features  then  led  to  the  selection  of  the  principal  parameters  and 
the  formulation  of  an  overall  design  procedure. 
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SECTION  2 

WIND  VELOCITY  DATA 


Although  atmospheric  soundings  have  been  made  for  many  years,  the  data 
usually  are  used  for  meteorological  purposes  for  which  they  were  originally 
intended.  It  was  not  until  1954  that  Norman  Slssenwine  (Ref.  7)  made  a  first 
attempt  at  defining  a  one  percent  wind  profile  for  use  in  missile  design 
(Fig.  3,  profile  No.  1).  This  profile  was  derived  from  sounding  data  avail¬ 
able  at  that  time,  with  particular  emphasis  on  the  North-Eastern  United  States, 
considered  the  windiest  part  of  the  continent.  The  approach  was  to  define, 
from  wind  statistics,  the  maximum  wind  velocity  (300  fps)  and  associated 
maximum  shear  (45  fps/ 1000  ft  »  0.  045  sec  *)  likely  to  be  exceeded  only  one 
percent  of  the  time  during  the  winter  season.  The  remainder  of  the  profile, 
above  and  below  the  shear  at  35,  000  ft  were  developed  using  ratios  of  the 
wind  velocities  of  the  limited  number  soundings  where  the  wind  velocity 
approached  300  fps  at  35,  000  ft  altitude. 

In  a  subsequent  study  (Ref.  8),  Slssenwine  developed  wind  profiles  for  Patrick 
AFB  at  various  probabilities  of  exceedance.  The  approach  is  very  nearly  the 
same  as  that  of  Ref.  7  with  the  exception  that  the  study  is  essentially  limited 
to  Patrick  AFB  sounding  data.  The  results  are  four  profiles  at  1,  5,  10,  and 
20  percent  calculated  risk  during  the  winter  time.  The  one  percent  profile 
(Fig.  3,  profile  No.  2)  has  a  maximum  wind  velocity  of  298  fps  at  45,  400  ft 
altitude,  accompanied  by  a  shear  of  0.  046  sec  *  and  shear  lengths  of  1000  ft 
below  and  3000  ft  above  the  peak.  The  remainer  of  the  profiles  consists  of 
actual  sounding  data  chosen  to  correspond  closely  with  the  maximum  wind 
velocities  and  shears. 

The  foregoing  studies,  as  well  as  one  conducted  by  Tolefson  (Ref.  2),  are 
all  based  on  data  collected  with  AN/GMD-1  sounding  equipment.  This 
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instrument  produces  considerable  inaccuracies  in  wind  velocity  measurements 
especially  at  high  altitudes  and  for  strong  wind  velocities!  resulting  in  gross 
inaccuracies  of  the  computed  wind  shears.  Detailed  error  studies  reported 
in  Refs.  1  and  2  show  that  the  RMS  error  is  often  of  the  same  order  of 
magnitude  as  typical  strong  shears.  In  Refs.  7  and  8,  Sissenwine  attempted 
to  account  for  these  errors  while  Tolefson  in  Ref.  2  did  not  correct  for  them 
in  the  shears  presented.  Consequently,  the  Sissenwine  method  tended  to 
reduce  extreme  shears  of  low  probability  and  Tolefson  tended  to  magnify 
these  values. 

With  the  development  of  the  AN/GMD-2  remitter  rawin  system,  instrument 
errors  were  reduced  by  about  one  order  of  magnitude  (Ref.  1).  Therefore, 
initial  measurements  with  this  system  reported  in  Ref.  1  are  extremely 
significant,  even  though  the  sample  is  small  by  statistical  standards.  It  was 
found  that  wind  shear  and  shear  length  are  related,  high  shears  being  associ¬ 
ated  with  thin  shear  layers.  For  example,  the  one  percent  maximum  shear 
is  0.  076  sec-*  for  1000  ft  shear  length  and  0.  050  sec  *  for  3000  ft  shear 
length.  An  important  feature  of  these  data  is  that  the  correlation  of  the 
probability  to  exceed  maximum  wind  velocity  and  shear  for  various  shear 
lengths  is  shown.  The  results  are  summarized  in  Table  1  and  plotted  in 
Fig.  2. 

Many  other  attempts  have  been  made  toward  the  definition  of  design  wind 
profiles.  Lockheed  devised  a  one  percent  profile  (Ref.  9)  consisting  of 
statistical  minimum  wind  velocity  with  a  transition  to  maximum  wind  velocity 
of  300  fps  in  the  altitude  interval  of  30,  000  to  40,  000  feet  (Fig.  4).  The 
transition  is  a  combination  of  a  0.  050  sec  *  wind  shear  over  2000  feet  and  a 
50  fps  gust  with  500  feet  gust  length.  Convair  (Ref.  11)  constructed  the  pro¬ 
file  shown  in  Fig.  4  and  Sissenwine  (Ref.  12)  revised  his  profile  as  illustrated 
in  Fig.  3  (profile  No.  3),  in  accordance  with  the  conclusions  reached  in 
Ref.  1. 

AviDyne  Research  conducted  extensive  studies  on  the  missile  response  to 
wind  shears  under  contract  AF  33(6l6)-5960  with  WADD  (Refs.  6  and  10). 
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A  statistical  sample  was  compiled  consisting  of  279  actual  wind  soundings 
considered  representative  of  the  seven  winters  from  which  it  was  derived. 

A  missile  was  "flown"  (on  the  computer)  through  each  one  of  these  profiles 
to  obtain  bending  moments  at  three  missile  stations.  Bending  moments  on 
the  missile  at  the  same  three  stations  were  obtained  by  subjecting  the  vehicle 
to  synthetic  profiles  and  statistical  arrays.  The  Lockheed  profile  appeared 
to  yield  the  best  results.  The  profile  recommended  for  design  in  Refs.  6  and 
10  Is  so  close  to  that  of  Lockheed  (Fig.  4)  that  it  is  not  plotted. 

Subsequent  studies  by  AviDyne  Research  (Ref.  13)  resulted  in  the  recommended 
profiles  shown  in  Fig.  5.  The  disadvantage  of  these  profiles  is  that  they  apply 
only  to  one  percent  probability,  leaving  the  designer  no  room  to  adjust  the 
profile  to  the  particular  mission  requirements  and  the  possible  limitations  of 
an  available  booster. 

A  summary  of  the  major  characteristics  of  the  wind  profiles  discussed  in  this 
chapter,  including  profiles  received  most  recently  from  Marshall  Space  Flight 
Center  (Fig.  6),  is  presented  in  Table  2. 
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SECTION  3 
VEHICLE  RESPONSE 


From  the  foregoing  it  is  apparent  that  the  selection  of  a  satisfactory  critical 
design  wind  profile  is  hard  to  accomplish  if  the  problem  is  approached  from 
the  sounding  data  only.  Reasonably  good  agreement  exists  on  the  maximum 
wind  velocity  and  shear  but  the  path  from  ground  to  shear  altitude  is  subject 
to  wide  speculation.  The  confusion  is  compounded  by  the  preferences  of  the 
various  monitoring  agencies:  The  Department  of  Defense  (Ref.  15)  specifies 
the  use  of  the  wind  profile  of  Ref.  8  (Fig.  3);  WADD  (Ref.  10)  recommends  a 
profile  nearly  identical  to  that  of  Ref.  9  (Fig.  4);  Aerospace  Corporation  has 
a  preference  for  the  profile  of  Ref.  13  (Fig.  5);  while  NASA  appears  to  have 
most  confidence  in  the  data  compiled  by  its  Marshall  Space  Flight  Center 
(Ref.  14,  Fig.  6). 

Actually,  the  maximum  response  of  a  vehicle  is  the  result  of  the  effects  of  the 
wind  profile  beneath  the  altitude  of  maximum  shear.  The  position  of  the 
vehicle  at  maximum  shear  is  dependent  upon  its  response  at  lower  altitudes. 
Therefore,  the  only  way  in  which  the  aerospace  industry  can  work  itself  out 
of  the  wind  profile  dilemma  appears  to  be  the  study  of  the  response  of  different 
vehicles  to  the  various  profiles. 

The  response  to  the  data  presented  in  Section  2  was  obtained  by  means  of  an 
IBM  7090  6-degree  of  freedom  trajectory  computer  program  which  incor¬ 
porated  the  characteristics  of  an  existing  missile  (Missile  A).  The  analysis 
is  limited  to  rigid  body  effects  since  winds  capable  of  being  measured  with 
present  radiosonde  equipment  do  not  appear  to  excite  flexible  body  modes 
(Ref.  6).  The  results  are  summarized  in  Table  3  and  typical  response  curves 
are  shown  in  Fig.  7.  The  analyses  were  conducted  for  three  wind  directions 
(head,  side,  and  tail),  disregarding  the  structural  and  performance 
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capabilities  of  the  missile.  In  some  head  wind  cases  this  necessitated  the  re¬ 
moval  of  the  engine  stops*  from  the  program.  Some  conclusions  are  obvious; 

a.  The  profiles  of  Refs.  9  and  12  (Lockheed  and  Sissenwine  respect¬ 
ively)  are  nearly  as  severe  with  the  Lockheed  profile  causing  the 
highest  response. 

b.  Although  ground  winds  vary  from  zero  to  35  fps,  their  effect  on 
the  maximum  response  is  not  immediately  apparent.  High  ground 
winds,  however,  appear  to  excite  a  long  period  mode  and  may  be 
useful  for  the  evaluation  of  damping  characteristics  of  the  vehicle. 

PARAMETER  STUDY 

To  isolate  the  principle  parameters  relating  wind  shear  and  vehicle  response, 
extensive  parameter  studies  were  conducted.  To  properly  describe  the 
response,  the  parameters  aq  and  0q  were  determined  to  be  definitive.  These 
parameters  correlate  closely  with  the  maximum  bending  moments  in  the 
missile  structure  and  since,  for  the  wind  directions  analyzed,  one  is  neg¬ 
ligible  when  the  other  is  critical,  the  correlation  is  with  resultant  bending 
moments.  It  evolves  that  the  response  to  each  wind  profile  is  defined  by 
one  value  of  these  parameters  regardless  whether  the  a  or  0  response  has 
one  maximum,  two  maxima,  or  low  long  period  damping. 

In  the  search  for  a  definitive  wind  profile  parameter,  numerous  variables 
were  investigated.  As  noted  above,  the  effect  of  ground  wind  is  not  obvious 
at  this  stage.  There  is  no  direct  relationship  between  wind  shear  and  vehi¬ 
cle  response.  Shear  length  affects  the  response  somewhat  but  not  in  a 
consistent  manner.  The  wind  velocity  increment  during  shear  is  related  to 
the  vehicle  response  but  the  effect  is  not  conclusive  enough.  Inconsistencies 
in  this  plot  all  had  in  common  that  this  variable  does  not  involve  the  effects 
of  the  wind  profile  beneath  the  altitude  of  maximum  shear.  At  this  point , 
logic  dictated  the  next  step:  The  integrated  effects  of  a  wind  profile  can  be 
described  by  the  maximum  wind  velocity  V  in  conjunction  with  the  inte¬ 
grated  area 

H 

f  cr  2 

A  =  J  VdH  ft  /sec 
o 

*  Restraints  built  into  the  engine  gimbal  mechanism  which  limit  the  engine 
deflections  to  a  predetermined  maximum  value. 
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under  the  profile,  where  Hcr  la  the  altitude  of  maximum  wind  velocity.  All 
profiles  with  the  same  maximuni  wind  velocity  (Vmax)  fall  on  a  single  curve 
of  maximum  response  (  a  q  or  0  q  )  versus  integrated  area  (A),  regardless 
of  the  wind  direction  (Fig.  8).  The  significance  of  the  area  A  is  underscored 
even  more  by  the  finding  that  a  change  in  critical  wind  direction  occurs  at  a 
specific  value  of  the  integrated  area.  For  vehicle  A  this  is  3.4  x  10^  ft^/sec, 
side  wind  being  critical  for  larger  values  of  A  and  head  wind  for  smaller 
values. 

In  order  to  Investigate  these  relationships  still  further,  a  variety  of  synthetic 
profiles  was  constructed  (Table  4,  Fig.  9)  and  the  vehicle  "flown"  through 
them  on  the  computer.  All  these  profiles  have  a  maximum  wind  velocity  of 
300  fps  and  were  constructed  around  the  one  percent  shears  presented  in 
Table  1: 

.075  sec'1  with  1000  ft  shear  length 
,  050  sec  1  with  3000  ft  shear  length 

It  is  seen  that  these  profiles  range  widely  up  to  excessively  severe.  The 
vehicle  response  for  these  profiles  is  plotted  in  Fig.  10.  Besides  the  obvious 
conclusion  that  the  data  trend  is  Identical  to  that  of  Fig.  8,  several  other 
observations  can  be  made: 

a.  Variations  of  critical  altitude  induce  response  changes  which 
are  parallel  with  or  on  the  conservative  side  of  the  general  data 
trend. 

b.  Increasing  ground  winds  cause  a  reduced  response  but  also  a 
shift  of  the  curve  into  the  critical  direction. 

The  possibility  of  establishing  a  family  of  curves  for  a  range  of  ground  wind 
probabilities  was  investigated.  However,  this  approach  was  abandoned  since 
the  present  Btate -of -th e-art  is  such  that  combined  probabilities  for  ground 
wind  and  high  altitude  shear  cannot  be  determined  with  any  degree  of  reliability 
(Ref.  Section  5).  In  the  following,  therefore,  ground  wind  (up  to  the  approxi¬ 
mate  one  percent  maximum  of  35  fps)  is  considered  as  a  random  variable 
which  contributes  to  the  data  scatter. 
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Further  investigation  showed  that  a  response  envelope  may  be  obtained  in 
parts,  one  for  each  wind  direction.  All  side  wind  data,  both  critical  and 
non- critical,  plot  as  an  approximate  straight  line  and  the  same  holds  for  all 
head  wind  data.  Consequently,  regression  analyses  were  performed  on  both 
sets  of  data  with  results  as  shown  in  Figs.  11  and  12.  The  high  value  of  the 
correlation  coefficient  (  r  =  0.  980)  is  mathematical  evidence  of  the  excellent 
linearity  of  the  data  (  r  =  1  for  perfect  linear  fit).  The  scatter  of  the  data 
about  the  mean  regression  line  is  expressed  by  the  standard  error  and  is 
suitably  small  (S^.  =  206  and  306  deg-psf  respectively).  Approximate  confi¬ 
dence  limits  at  one -standard  error -true  confidence  limits  are  hyperbolic  - 
contains  the  data  with  only  a  few  exceptions.  A  composite  of  the  head  and 
side  wind  regressions  is  presented  in  Fig.  13.  In  order  to  obtain  an  enve¬ 
lope  of  the  data,  the  confidence  interval  for  the  head  wind  data  is  taken  as 
1.41  standard  errors.  It  is  seen  that  only  the  response  of  profile  S21  exceeds 
the  envelope.  The  probability  of  conditions  represented  by  this  profile  ever 

occurring  (Fig.  9)  is  so  remote  that  it  will  be  disregarded.  The  linear 

6  2 

components  of  the  response  envelope  intercept  at  A  =  3.  385  x  10  ft  /sec 
which  is  in  excellent  agreement  with  a  previous  observation  based  upon  the 
individual  data. 


DESIGN  PROFILES 


The  foregoing  procedure  for  obtaining  a  response  envelope  is  entirely  too 
complex  and  laborious  to  be  used  in  the  design  and  launch  control  of  missile 
systems.  Instead,  the  statistical  envelope  should  be  approximated  by  a 
minimum  number  of  simple  profiles.  Inspection  of  Figs.  11,  12,  and  13 
Indicates  that  out  of  all  the  wind  profiles  studied,  there  are  four  that  are 
definitive  (See  also  Section  V).  These  are  (Ref.  Tables  2  and  4): 

a.  The  revised  Sissenwine  profile  (No.  3)  which  is  based  entirely 
upon  the  one  percent  AN/GMD-2  data  of  Ref.  1  (Table  1).  This 
profile  is  constructed  with  1000,  3000,  and  5000  ft  shears 
(0.  075,  0.  050,  and  0.  033  sec-1  respectively)  that  are  associated 
with  a  maximum  wind  velocity  of  300  fps.  The  critical  altitude 
is  30,000  ft  (Fig.  3). 
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b.  Profile  No.  3c  which  is  the  same  as  the  revised  Sissenwine 
profile  (No.  3)  except  for  a  critical  altitude  of  40,000  ft. 

c.  The  minimum  area  profile  Dl,  300.  This  high-intensity  shear 
profile  has  an  integrated  area  that  may  well  be  smaller  than  is 
physically  attainable  and,  therefore,  results  in  a  response  that 
can  hardly  be  exceeded.  The  area  was  minimized  by  the  selection 
of  a  low  critical  altitude  (about  30,000  ft)  and  extension  of  the 
1000  ft  shear  associated  with  the  maximum  wind  velocity  — 

0.  075  sec-1  for  y  =  300  fps  (Table  1)  —to  3000  ft  shear  length 
(Fig.  9).  maX 

d.  The  maximum  area  profile  D2,  300.  The  area  of  this  profile  was 
maximized  by  selecting  a  high  critical  altitude  (about  40,  000  ft) 
and  decreasing  the  wind  velocity  linearly  from  maximum  to  zero 
ground  wind.  (Fig. 9. ) 

The  first  two  profiles  (Nos.  3  and  3c)  form  a  pair  that  is  relatively  simple 
and  is  rationally  derived  from  recent  wind  shear  data  (Table  1).  However, 
the  area-range  of  this  pair  is  relatively  small  and  is  exceeded  by  several 
other  literature  wind  profiles  (Fig.  8).  In  the  following  they  will  be  treated 
as  design  profiles  and  denoted  by  D3,  300  and  D4,  300  respectively. 

The  extreme  area  profiles  Dl,  300  and  D2,  300  are  related  to  the  statistical 
wind  data  Table  1.  They  bracket  the  range  of  integrated  areas  that  is 
physically  attainable  for  a  specific  maximum  wind  velocity.  Since  both 
profiles  are  already  extreme,  the  effect  of  a  finite  ground  wind  has  been 
purposely  ignored. 

Since  all  four  design  profiles  are  related  to  the  same  statistical  wind  shear 
data  (Table  1),  they  may  readily  be  expanded  into  series  (i.  e.  Dl,300; 

Dl,  250;  Dl,  225;  etc.  )  with  maximum  wind  velocity  as  the  argument.  For 
instance,  design  profile  D3,  225  is  constructed  similar  to  profile  D3.300 
except  that  the  shears  associated  with  the  maximum  wind  velocity  of  225  fps 
(Table  1)  are  used.  The  detail  development  of  each  one  of  the  four  series 
is  presented  in  Appendix  A.  In  order  to  facilitate  the  direct  utilization  of 
any  profile  in  the  series,  the  coordinates  and  integrated  areas  are  presented 
in  tabular  and  graphical  form. 
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EFFECT  OF  DIFFERENT  VEHICLES 


The  development  of  a  design  procedure  would  be  premature  at  this  point 
since  all  the  data  considered  so  far  were  associated  with  only  one  vehicle 
(A).  In  order  to  obtain  the  broadest  possible  base  for  the  design  procedure^ 
two  other  vehicles  (B  and  C)  were  selected  for  analysis.  The  three  vehicles 
A,  B,  and  C  are  as  radically  different  as  is  feasible  at  present.  The  two 
additional  vehicles  were  "flown"  through  the  four  design  profile  series,  as 
well  as  a  selected  number  of  the  individual  wind  profiles  used  before  on 
vehicle  A.  The  results  for  all  three  missiles  systems  are  shown  in  Figs. 

14,  15,  and  16.  It  may  be  seen  that  the  same  trends  hold  for  all  three  vehi¬ 
cles  and  design  profile  series  D1  through  D4  always  define  the  response 
envelope.  This  is  so  because  they  approximate  the  extremes  and  the  mode 
of  the  statistical  distribution  of  integrated  area.  Note  also  that  a  maximum 
wind  velocity  of  about  250  fps  appears  to  be  optimum;  greater  maximum 
wind  velocities  may  be  very  costly  in  terms  of  structural  weight  without 
commensurate  gains  in  probability. 

SUMMARY  OF  PARAMETRIC  RELATIONSHIPS 


Prior  to  the  development  of  a  design  procedure,  a  summary  of  the  various 
parameters  and  their  relationships  is  warranted: 

a.  The  basic  parameters  that  define  the  response  to  winds  aloft  of 
missile  systems  studied  are: 

(1)  Maximum  wind  velocity) 

(2)  Integrated  area  Jo  CT  Vw  dH  where  Hcr  is  the  altitude  at 
which  the  maximum  wind  velocity  first  occurs. 

This  is  valid  regardless  of  the  missile  system  and  the  severity 
of  the  wind  shears. 

b.  The  maximum  response  for  all  conceivable  wind  profiles  having 
the  same  maximum  wind  velocity  is  a  well-defined  function  of 
the  integrated  area.  This  relationship  holds  not  only  when  the 
wind  direction  is  held  constant  (Figs.  11  and  12),  but  also  when 
only  critical  wind  directions  are  considered  (Fig.  13). 
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c.  Within  statistical  limits,  the  vehicle  response  decreases  with 
increasing  integrated  area. 

d.  Response  envelopes  (constant  maximum  wind  velocity,  critical  wind 
directions,  Figs.  13  to  16)  exceed  the  responses  associated  with  all 
wind  profiles  (whether  obtained  from  the  literature  (Table  2)  or 
arbitrarily  (Table  4))  by  no  more  than  4.5  percent  on  the  average 
within  a  range  from  0  to  10  percent. 

e.  For  constant  maximum  wind  velocity  there  are  two  critical  wind 
directions,  side  wind  and  head  wind.  Of  these,  side  wind  is  by 
far  the  predominant  one.  However,  it  is  conceivable  that  for 
certain  pitch  programs  tall  wind  might  be  a  critical  direction 
rather  than  head  wind. 

f.  Design  profiles  D1  through  D4  define  the  response  envelope  for  a 
specific  maximum  wind  velocity.  Because  they  are  based  upon 
statistical  wind  data  (Table  1),  they  are  expandable  into  series 
with  maximum  wind  velocity  as  the  argument  (Appendix  A). 

g.  Design  profiles  D1  and  D2  are  extreme  area  profiles  which  may 
be  considered  as  boundary  conditions  for  the  area-response  func¬ 
tion.  They  result  in  wedge-type  boundaries  within  which  the 
response  envelopes  are  contained  (Figs.  14,  15,  and  16). 

h.  Design  profiles  D3  and  D4  are  made  up  of  statistical  wind  shears 
based  on  actual  sounding  data  (Ref.  1).  They  should  result, 
therefore,  in  a  probably  or  most  likely  response.  This  is  veri¬ 
fied  by  the  results;  the  response  to  profiles  D3  and  D4  (or  3  and 
3c)  is  centered  between  the  limits  whenever  they  are  critical 
(Figs.  11  through  16). 

i.  Design  profile  D2  is  more  realistic  than  might  appear  at  first 
glance:  Its  correlation  throughout  with  Sissenwine  profile  No.  2 
is  striking  (Figs.  14  to  16). 
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SECTION  4 
DESIGN  PROCEDURE 

Within  the  limits  of  statistical  confidence,  the  response  of  a  missile  system 
to  winds  aloft  is  completely  defined  by  the  diagrams  of  Figures  14  to  16.  The 
diagrams  include  the  response  to  winds  only  and  the  effects  of  elastic  response 
were  not  considered  based  on  conclusions  reached  in  Ref.  6.  For  design  pur¬ 
poses,  these  diagrams  are  rather  complicated  and  the  necessity  of  always 
having  to  know  the  integrated  area  makes  them  laborious.  Therefore,  a 
simplified  approach  would  be  well  worth-while. 

DEVELOPMENT 

One  of  the  conclusions  derived  from  the  response  diagrams  of  Figures  14  to  16 
(Section  3,  Page  12)  is  that  the  possible  response  is  bracketed  by  design  pro¬ 
files  D1  and  D2  while  the  probable  response  is  associated  with  design  — 
profiles  D3  and  D4  (see  also  Section  5).  From  this  it  follows  that  the 
integrated-area-parameter  can  be  eliminated  by  plotting  the  response  as  a 
function  of  maximum  wind  velocity  along  the  contour  of  a  design  profile  series. 
This  results  in  a  wedge-shaped  "scatter  band"  with  a  maximum  wind  velocity 
as  the  independent  variable  (Figs.  17,  18,  and  19).  The  band  is  bounded  by 
the  response  associated  with  design  profiles  D1  and  D2  while  the  most  likely 
response  associated  with  design  profiles  D3  and  D4  is  centrally  located  as 
expected.  The  significance  of  this  diagram  lies  in  the  following: 

a.  The  independent  variable,  maximum  wind  velocity,  is  a  direct 
variable  (directly  associated  with  the  physical  data)  rather  than 
a  derived  variable  like  the  integrated  area. 

b.  The  diagram  may  be  directly  associated  with  probability  of 
occurrence  by  visualizing  a  probability  scale  in  the  third  dimen¬ 
sion.  The  bell- shaped  curve  would  have  its  mode  between  D3  and 
D4  and  reduce  to  near- zero  at  D1  and  D2  (Ref.  Section  5). 
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The  design  response  curves  of  Figs.  17  to  19  do  not  negate  the  significance  of 
the  integrated  area  as  a  definitive  parameter.  To  the  contrary,  only  by  virtue 
of  the  integrated-area-parameter  can  it  be  stated  that,  for  a  specific  maximum 
wind  velocity  and  regardless  of  the  wind  profile,  the  missile  system  response 
should  not  exceed  D1  and  will  occur  with  the  greatest  likelihood  between  D3 
and  D4  (Ref.  Section  5). 

The  final  design  diagrams  are  presented  in  Figs.  20,  21  and  22.  They  con¬ 
sist  of: 

a.  The  design  wind  response  of  Figs.  17  to  19.  In  most  cases  there 
will  be  no  need  to  include  the  lower  boundary  (D2). 

b.  The  probability  to  exceed  the  maximum  wind  velocity  as  obtained 
from  Table  1.  Note  that  the  probability  curve  is  associated  only 
with  the  most  likely  response  (D3  and  D4),  never  with  the  envelope 
(Dl).  This  derives  from  the  observations  made  about  probability  of 
occurrence  in  (b)  above  (See  also  Section  5,  Page  22). 

c.  Any  capability  limitations  of  the  missile  system  such  as  structural, 
engine  gimbal  stops,  etc. 

The  design  diagram  may  be  put  to  many  uses  as  shown.  For  a  completed  mis¬ 
sile  system,  the  launch  probability  and  the  safe  maximum  wind  velocity  (with¬ 
out  detailed  analysis  of  pre-launch  wind  soundings)  may  be  determined  (Fig.  20). 
In  the  early  stages  of  design,  required  limit  capabilities  may  be  obtained  based 
upon  launch  probabilities  desired  by  the  customer  (Fig.  21).  Note  that  decreas¬ 
ing  the  probability  to  launch  from  say  99%  to  95%  may  result  in  relatively  vast 
weight  savings  and  that  this  relationship  is  nonlinear.  For  operational  missile 
systems  the  diagram's  usefulness  is  that  it  facilitates  a  pre-launch  check 
whether  limit  capability  may  be  exceeded  (Fig.  22).  A  detailed  trajectory  analy¬ 
sis  to  ascertain  of  safe  launch  conditions  is  required  only  when  the  maximum 
wind  velocity  obtained  from  wind  soundings  is  in  the  grey  area  where  limit 
capability  lies  between  the  most  likely  and  extreme  response. 

DETAILED  PROCEDURE 

The  actual  procedures  to  be  used  for  a  new  design  will  vary  with  the  require¬ 
ments  peculiar  to  the  missile  system.  Although  a  more  or  less  complete 
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procedure  is  outlined,  a  much  more  simplified  version  will  be  adequate  in 
most  instances.  For  instance,  in  some  cases  the  design  response  envelope 
(Dl)  may  be  all  that  is  required.  In  other  cases  only  one  probable  response 
curve  may  be  desired  (D3  or  D4,  whichever  is  maximum).  In  order  to  estab¬ 
lish  all  the  curves  of  the  design  diagram,  the  following  procedure  is  most 
efficient: 


a.  If  applicable,  update  Table  1  with  the  latest  data  from  the  literature 
and  revise  the  tables  in  Appendix  A  accordingly. 

b.  Perform  a  trajectory  analysis  for  the  no-wind  condition.  Study  the 
results  and  make  any  changes  required  in  the  control  and  command 
data. 

c.  Establish  the  extreme  maximum  wind  velocity  desired  (300  fps  in 
this  report),  construct  design  profile  Dl  and  perform  trajectory 
analyses  using  all  wind  directions  and  with  all  capability  limitations 
(i.  e.,  engine  stops)  removed.  Determine  the  two  most  critical  wind 
directions  and  perform  all  future  trajectory  analyses  with  these  two 
only. 

d.  Construct  design  profile  Dl  for  two  or  three  intermediate  maximum 
wind  velocities  and  perform  trajectory  analyses.  Start  the  design 
diagram  with  the  response  envelope  and  probability  curve  (from 
Table  1).  In  some  cases  this  may  be  all  that  is  required. 

e.  Determine  whether  the  probable  response  is  desired  for  both  design 
profiles  D3  and  D4  or  only  the  maximum  of  the  two.  Determine  the 
critical  altitude  as  required,  using  the  results  from  previous  anal¬ 
yses)  if  necessary,  perform  trajectory  analyses  for  the  extreme 
maximum  wind  velocity  at  various  critical  altitudes. 

f.  Construct  design  profiles  D3  and/or  D4  for  the  extreme  maximum 
wind  velocity  and  two  or  three  intermediate  values  and  perform 
trajectory  analyses.  Plot  the  probable  response  on  the  design 
diagram. 

g.  If  desired,  construct  design  profile  D2  using  no  more  than  three 
maximum  wind  velocities  Including  the  extreme  one. 

h.  Complete  the  design  diagram  with  all  capability  limitations  that  may 
prevail.  This  includes  any  limitations  that  were  temporarily  re¬ 
moved  from  the  trajectory  analysis  program. 


EFFECT  OF  WIND  DIRECTION 


At  times,  a  missile  system  is  specifically  intended  for  launching  from  one 
base  only.  If  this  is  a  base  with  a  predominant  wind  direction  such  as 
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Cape  Canaveral,  weight  savings  may  be  effected  by  accounting  for  this  bias. 
Rather  than  to  assume  the  maximum  wind  velocity  to  be  constant  regardless 
of  the  wind  direction  as  was  the  basis  for  the  design  procedure  presented 
above,  the  maximum  wind  velocity  may  be  varied  with  the  wind  direction  in 
accordance  with  local  wind  statistics.  One  way  to  accomplish  this  is  by 
approximating  the  wind  velocity  at  a  certain  altitude  and  for  the  desired  prob¬ 
ability  of  occurrence  by  an  ellipse  with  axes  in  the  North-South  and  East-West 
directions.  With  the  launch  azimuth  known,  the  maximum  head.  Bide,  and  tail 
wind  may  then  be  scaled  off  the  ellipse.  For  Cape  Canaveral,  for  instance, 
this  procedure  results  in  approximate  maximum  wind  velocities  as  follows: 

151  fps  head  wind 
183  fps  side  wind 
300  fps  tail  wind 

Considering  that  head  wind  and  side  wind  usually  are  critical,  this  procedure 
may  result  in  considerable  weight  savings.  Extreme  caution  must  be  exer¬ 
cised,  however,  since  a  missile  system  designed  to  these  limited  criteria 
may  not  be  structurally  adequate  for  any  other  launch  site  than  the  one  used 
as  a  basis  for  design. 

During  pre-launch  wind  soundings,  it  has  sometimes  happened  that  the  maxi¬ 
mum  wind  velocity  at  altitude  was  moderate  but  its  direction  changed  quite 
drastically.  This  is  probably  due  to  two  layers  of  air  moving  into  different 
directions.  In  such  a  case,  the  missile  system  response  should  be  greater 
than  that  estimated  from  the  measured  maximum  wind  velocity  using  the 
design  diagram  (Fig.  22).  How  to  account  for  this  possibility  should  be  the 
subject  of  a  future  investigation. 

PRE  -  LAUNCH  CHECKOUT  PROCEDURE 

The  procedure  developed  in  this  report  is  eminently  suited  to  facilitate  pre¬ 
launch  checkout.  It  can  be  used  at  the  launch  site  rather  than  at  some  far¬ 
away  control  center.  The  only  requirement  is  that  some  facility  be  available 
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for  the  progressive  Integration  of  the  wind  velocity  data  being  obtained  from 
radiosonde  measurements.  The  detailed  procedure  would  be  as  follows: 

a.  For  the  particular  vehicle  being  launched,  two  plots  should  be  avail¬ 
able.  The  first  one  should  be  of  the  type  shown  in  Figs.  14  to  16  and 
the  other  one  like  Figs.  20  to  22. 

b.  From  radiosonde  measurements,  the  maximum  wind  velocity  and 
associated  integrated  area  are  obtained.  The  wind  direction  should 
be  checked  to  ensure  that  there  is  no  drastic  change  resulting  in 
two  high  shears. 

c.  On  the  first  plot  (like  Figs.  14  to  16),  a  vertical  line  is  drawn  at  the 
Integrated  area  associated  with  the  measured  maximum  wind  velocity. 
Then  an  envelope  is  estimated  for  the  measured  maximum  wind  velo¬ 
city.  The  intersect  of  integrated  area  and  maximum  wind  velocity 
envelope  is  then  located  with  respect  to  the  design  profile  series,  just 
above  D4  for  instance.  Note  that  the  response  is  not  read  from  this 
plot  since  the  maximum  wind  velocity  envelope  is  only  estimated  and, 
therefore,  inaccuracies  would  be  compounded. 

d.  Next,  the  second  plot  (like  Figs.  20  to  22)  is  entered  with  the  meas¬ 
ured  maximum  wind  velocity  and  a  vertical  line  drawn  into  the 
response  band  up  to  the  location  obtained  under  (c  ),  just  above  D4 
in  this  case.  The  response  associated  with  this  point  may  now  be 
read  and  structural  and  control  margins  computed. 

e.  The  last  step  should  be  to  estimate  the  launch-probability  associated 
with  the  radiosonde  data  and  to  compare  it  with  the  design  probability 
level.  The  difference  between  the  two  is  an  additional  measure  for 
the  degree  of  safety  associated  with  the  particular  launch. 


The  foregoing  procedure  has  actually  been  used  during  several  recent  launches 
and  was  cross-checked  against  6D  trajectory  computer  runs  for  the  overall 
wind  profile  as  measured  by  radiosonde.  The  correlation  to  date  has  been 
excellent. 

Further  improvements  in  this  procedure  may  be  obtained  in  the  future  when 
statistics  of  the  integrated  area  become  available.  Then  the  response  band 
will  be  defined  by  curves  for  specific  probabilities  to  exceed,  rather  than 
design  profile  series  D1  through  D4.  As  a  result,  the  joint  probability  of 
measured  maximum  wind  velocity  and  associated  Integrated  area  may  be 
estimated  and  the  response  expressed  in  terms  of  probability  to  exceed.  This 
will  provide  a  more  rational  basis  for  the  prediction  of  the  expected  response 
during  the  actual  launch. 
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SECTION  5 

PROBABILITY  CONSIDERATIONS 

In  the  foregoing!  relatively  little  was  said  about  probabilities  of  occurrence 
and  in  the  development  of  the  design  diagram  the  probability  concept  was 
without  further  discussion  associated  with  maximum  wind  velocity  only.  This 
was  done  for  purposes  of  clarity  and  in  the  following  the  various  statistical 
aspects  will  be  discussed  in  broad  terms  only. 

WIND  PROFILES 

The  general  approach  to  date  has  been  to  construct  a  synthetic  wind  profile 
for  a  certain  probability  of  occurrence  using  statistical  wind  sounding  data. 
The  result  is  a  large  variety  of  wind  profiles  (Table  2)  which  includes  re¬ 
sponses  that  vary  widely  (Fig.  8).  The  major  cause  for  this  wide  spread  is 
that  the  construction  of  a  synthetic  wind  profile  for  a  specific  probability  of 
occurrence  is  highly  subjective  and  very  nearly  impossible  to  accomplish 
rigorously.  Even  if  it  is  assumed  that  all  the  parameters  that  define  a  wind 
profile  individually  are  accurately  defined  in  terms  of  probability  of  occur¬ 
rence  -  which  they  are  not  -  then  the  probability  level  of  the  composite  profile 
is  still  a  matter  of  conjecture.  An  example  will  illustrate  this. 

The  requirement  is  to  define  a  wind  profile  with  a  probability  to  be  exceeded 
of  one  percent.  Let  it  be  assumed  that  the  AN/GMD-2  data  of  Table  1  are 
statistically  reliable  estimates.  Further  assume  that  the  one  percent  ground 
wind  is  35  fps.  The  composite  constitutes  the  revised  Sissenwine  profile 
No.  3  (Table  2)  or  design  profile  D3,  300.  The  average  probability  for  this 
composite  to  be  exceeded  is  composed  of  probability  functions  associated  with 
the  maximum  wind  velocity,  1000  ft  shear,  3000  ft  shear,  5000  ft  shear, 
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35  fps  gound  wind,  and  critical  altitude.  While  maximum  wind  velocity  and 
shears  individually  are  one  percent  values,  the  probability  of  these  occurring 
simultaneously  is  at  best  60  percent  (Ref.  1),  resulting  in  a  joint  probability 
less  than  one  percent.  All  evidence  indicates  that  winds  aloft  and  ground 
winds  are  independent  variates  (statistical  variable),  meaning  that  their 
probability  functions  must  be  multiplied  to  obtain  their  joint  probability 
function.  The  result  is  an  effective  probability  to  exceed  vastly  less  than 
one  percent.  Last  but  not  least,  the  critical  altitude  at  which  the  maximum 
wind  velocity  occurs  has  its  own  independent  distribution.  The  probability 
for  the  critical  altitude  to  occur  specifically  at  the  altitude  of  maximum 
response  of  a  certain  missile  system  is  always  very  small  (cross-hatched 
area  in  sketch).  In  summary,  the  probability  associated  with  a  composite 
wind  profile  can  be  defined  only  in  terms  of  "much  less  than,  "  much  less 
than  one  percent  in  the  above  example. 


DESIGN  DIAGRAM 

Even  if  highly  reliable  statistical  data  were  available  on  all  the  parameters 
that  enter  into  the  vehicle  response  to  winds  aloft,  a  project  to  compute  the 
associated  probabilities  would  not  be  worth-while.  The  reason  is  that  the 
result  would  be  average  probabilities  to  exceed  which  in  turn  have  their  own 
statistical  distributions.  No  matter  how  great  the  computational  effort,  the 
response  for  a  certain  probability  level  may  only  be  determined  within 
confidence  limits.  The  effect  of  increasing  analysis  effort  would  only  be  to 
narrow  these  confidence  limits.  Since  the  physical  phenomenon  is  a  random 
process,  rather  wide  confidence  limits  are  inherent  of  the  data. 
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With  the  design  procedure  developed  in  Section  4,  an  attempt  has  been  made 
to  account  for  this  difficulty  of  statistical  definition.  The  discussion  in  the 
preceding  paragraphs  may  also  be  reversed  and  the  statement  made  that  all 
the  wind  profiles  presented  in  the  literature  (Section  2)  are  correct  in  that 
each  one  of  them  constitutes  a  portion  of  the  overall  phenomenon.  While 
most  of  these  profiles  are  in  agreement  with  regard  to  the  one  percent  and 
five  percent  maximum  wind  velocities,  their  shears  vary  widely  as  is  shown 
conclusively  in  Figs.  23  and  24.  Nonetheless,  their  response  as  a  function 
of  integrated  area  is  consistent  within  rather  narrow  limits  (Figs.  8  and  13). 
This  leads  to  the  inevitable  conclusion  that  the  maximum  wind  velocity  is  a 
principal  parameter  and  the  effect  of  the  other  variables  (not  including  inte¬ 
grated  area)  is  limited  to  a  contribution  to  the  dispersion. 

In  studying  plots  like  Figs.  13  and  14,  two  types  of  dispersion  are  of  concern: 
(a)  of  the  response  about  the  mean  regression  curve  and  (b)  the  wide  spread 
of  the  integrated  area. 

a.  Deviations  from  the  mean  regression  curve  are  caused  by  secondary 
variables  such  as  critical  altitude,  ground  wind,  and  actual  shears. 
The  effect  of  critical  altitude  and  ground  wind  have  already  been 
discussed  (Fig.  10).  A  study  of  individual  data  points  (Table  2, 

Fig.  8)  indicates  that  for  constant  integrated  area  low  average  shears 
cause  a  response  on  the  lower  side  of  the  confidence  band  (profile 
nos.  5a,  5b).  However,  since  the  maximum  deviation  from  the 
envelope  is  only  of  the  order  of  10  percent,  refined  techniques  for 
the  evaluation  of  the  secondary  variables  do  not  appear  to  be 
worthwhile. 

b.  Next  to  maximum  wind  velocity,  the  integrated  area  is  the  basic 
parameter  which  defines  the  vehicle  response.  Its  range  derives 
from  all' the  different  wind  profiles  which  are  possible  at  a  constant 
maximum  wind  velocity  (Fig.  9).  These  profiles  vary  randomly, 
some  of  them  occurring  more  frequently  than  others.  Therefore, 
the  Integrated  area  has  a  statistical  distribution  which  is  peculiar  to 
the  physical  phenomenon  under  investigation.  Its  distribution  may 
be  determined  directly  from  the  wind  sounding  data  by  simple  data 
reduction  methods.  Because  of  its  pronounced  effect  upon  the 
vehicle  response,  a  detailed  knowledge  of  the  statistical  distribution 
of  this  parameter  is  highly  desirable. 

Since  the  statistical  distribution  of  the  integrated  area  is  not  known  at  present, 
its  effect  must  be  assessed  qualitatively.  Referring  to  Figs.  14  to  16,  the 
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probability  distribution  of  the  integrated  area  may  be  visualized  in  the  third 
dimension.  Along  the  300  fps  envelope,  the  probability  is  near  zero  at  Dl, 
rises  to  a  maximum  somewhere  around  D3  and  D4,  and  drops  off  to  near  zero 
at  D2.  Similar  observations  can  be  made  in  Figs.  17  to  19. 

In  summary,  the  selection  of  design  profiles  D3  and  D4  to  induce  the  most 
likely  vehicle  response  is  based  upon  the  following  considerations: 

a.  The  statistical  distribution  of  the  integrated  area  should  have  a 
mode  somewhere  around  03  and  D4. 

b.  Design  profiles  03  and  D4  are  based  upon  the  most  complete 
statistical  wind  shears  published  to  date  (Ref.  1). 

c.  The  basic  design  profiles  D3,  300  (No.  3)  and  D4,  300  (No.  3c)  are 
expandable  into  series  which  permits  the  determination  of  design 
profiles  for  any  desired  maximum  wind  velocity  (Appendix  A). 

Once  the  distribution  of  the  integrated  area  is  known,  it  is  possible  that  a 
simplified  procedure  may  make  the  use  of  the  design  profiles  superfluous. 

Until  such  time,  however,  the  design  profiles  constitute  the  only  means  for 
estimating  the  mode  and  limits  of  the  response  variation  due  to  integrated 
area. 

The  reasons  for  the  format  of  the  design  diagram  (Figs.  20  to  22)  should  now 
be  clear.  Reasonable  reliable  statistics  are  available  for  the  maximum  wind 
velocity  and,  therefore,  are  shown.  No  statistics  are  available  for  the  inte¬ 
grated  area  and  its  effect  is  estimated  with  a  response  envelope  (design 
profiles  Dl)  and  a  probable  response  band  (design  profiles  03  and  D4). 

In  using  the  design  diagrams,  it  should  be  borne  in  mind  that  the  probability 
to  exceed  the  lower  envelope  (D2)  is  about  100- percent,  the  probability  to 
exceed  the  mean  response  (around  03  and  04)  is  50  percent,  and  the  probabil¬ 
ity  to  exceed  the  upper  envelope  is  close  to  zero.  Therefore,  the  joint 
probabilities  at  a  one  percent  maximum  wind  velocity  for  instance  are  approx¬ 
imately  one  percent,  0.5  percent,  and  zero  respectively.  It  follows  that  the 
recommended  design  procedure  to  relate  probability  with  the  probable  response 
band  is  conservative. 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


An  elaborate  investigation  of  high  altitude  winds  and  their  effect  upon  verti¬ 
cally  rising  missile  systems  resulted  in  the  isolation  of  two  definitive  param¬ 
eters,  maximum  wind  velocity  and  integrated  area  under  the  wind  profile. 
Within  statistical  limits,  these  two  parameters  completely  define  the  vehicle 
response.  For  constant  maximum  wind  velocity,  the  response  is  a  monoto¬ 
nously  decreasing  function  of  integrated  area.  All  "design"  wind  profiles 
recommended  in  the  literature  are  discrete  cases  of  this  general  relationship. 
For  the  class  of  missile  systems  Btudied,  the  vehicle  response  to  winds  aloft 
needs  to  be  determined  for  a  rigid  body  only. 

The  probability  of  a  certain  response  to  be  exceeded  depends  jointly  upon  the 
probability  functions  of  maximum  wind  velocity  and  integrated  area.  For  lack 
of  statistics  for  the  integrated  area,  this  joint  probability  at  present  must  be 
estimated.  This  is  accomplished  with  the  aid  of  design  profiles  described  in 
detail  in  Appendix  A.  The  specific  function  of  these  design  profiles  is  as 
follows:  Design  profiles  D1  and  D2  define  response  envelopes  associated 
with  probabilities  to  exceed  integrated  are?,  of  approximately  zero  percent 
and  100  percent  respectively.  Design  profiles  D3  and  D4  define  the  probable 
response  band  associated  with  an  estimated  mean  of  the  integrated  area 
distribution.  Theoretically,  the  probability  to  exceed  a  certain  response  is 
the  product  of  the  distribution  functions  of  maximum  wind  velocity  and 
integrated  area.  This  means  that  for  a  one  percent  maximum  wind  velocity, 
design  profile  D2  (lower  envelope)  induces  a  response  that  may  be  exceeded 
about  one  percent  of  the  time,  design  profiles  D3  and  D4  (probable  response) 
approximately  0.5  percent  of  the  time,  and  design  profile  D1  (upper  envelope) 
close  to  zero  percent  of  the  time.  As  long  as  adequate  statistical  data  are  not 
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available,  however,  the  conservative  approach  is  recommended  to  associate 
maximum  wind  velocity  probabilities  with  the  probable  response.  In  the  above 
example  this  means  that  design  profiles  D3  and  D4  are  assumed  to  induce  a 
response  that  may  be  exceeded  one  percent  of  the  time. 

The  results  of  this  investigation  were  synthesized  in  a  design  diagram  (Figs. 
20,  21,  and  22)  which  is  based  upon  the  parametric  relationships  and  includes 
the  probability  considerations  discussed  above.  A  detailed  procedure  for  the 
construction  of  this  diagram  has  been  presented.  The  design  diagram  may  be 
used  for  a  variety  of  purposes  during  all  phases  of  the  design,  including  a 
check  on  the  missile  system's  capability  to  withstand  winds  measured  prior 
to  launch. 


A  major  advantage  of  the  proposed  procedure  is  its  flexibility  which  permits 
the  designer  to  trade-off  launch-probability  against  cost  in  weight  and  dollars. 
Unlike  with  most  existing  wind  criteria,  it  is  now  possible  to  back-off  from  an 
initial  launch-probability  which  design  analysis  may  prove  1o  be  too  high,  while 
retaining  a  reasonably  accurate  knowledge  of  performance  gains  and  the  re¬ 
sulting  launch  probability. 


Some  miscellaneous  conclusions  of  interest  are: 

a.  The  relationship  between  vehicle  response  and  probability  of 
occurrence  is  nonlinear.  As  a  consequence,  increasing  the 
probability  to  launch  by  a  few  percentage  points,  say  from  95  to 
99  percent,  will  cause  the  weight  to  increase  disproportionately. 

b.  If  a  missile  system  is  being  designed  for  launch  from  one  base  only 
and  if  this  is  a  base  with  a  predominant  wind  direction,  then  signifi¬ 
cant  weight  savings  may  be  accomplished  by  accounting  for  this  bias 
in  the  construction  of  the  design  diagram. 

c.  High  ground  winds  may  excite  a  long  period  mode  of  the  missile 
system  and  thus  are  useful  for  the  evaluation  of  damping  character¬ 
istics  of  the  vehicle. 

This  investigation  has  resulted  in  a  systematic  and  parametric  approach  for 
the  determination  of  loads  due  to  winds  aloft.  However,  the  overall  procedure 
is  incomplete,  principally  for  lack  of  data  and,  therefore,  has  purposely  been 
kept  conservative.  In  order  to  make  it  possible  to  obtain  more  accurate  loads 
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and  further  weight  reduction,  future  studies  in  the  following  areas  are  urgently 
recommended: 

a.  The  available  wind  sounding  data  should  be  re-evaluated  to  obtain 
statistics  of  the  integrated  area  under  the  wind  profile  up  to 
maximum  wind  velocity  at  critical  altitude.  These  statistics  should 
be  determined  as  a  function  of  maximum  wind  velocity.  It  would  be 
advantageous  if  the  statistical  distribution  of  critical  altitude  were 
to  be  obtained  simultaneously. 

b.  The  effect  of  drastic  changes  in  wind  direction  upon  the  vehicle 
response  should  be  studied  and  the  results  incorporated  in  the 
design  diagram. 


For  vertically  rising  missile  systems,  loads  due  to  gust  generally  are  less 
severe  than  wind-induced  loads.  However,  they  may  still  be  appreciable  and 
should  be  considered.  Since  any  gust-load  analysis  must  include  flexible  body 
effects,  its  consideration  is  beyond  the  scope  of  this  report.  Instead,  it  will 
be  made  the  subject  of  a  future  investigation. 
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Table  3 

MAXIMUM  RESPONSE  FOR  VARIOUS  WIND  PROFILES*  VEHICLE  A 


PN°o^ 

Fig. 

No. 

Tail  Wind 

Head  Wind 
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e  Wind 

deg 
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deg  -  psf 
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deg 
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deg 
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4.  89 
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6.  71 
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4.  82 
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4.  39 
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5.03 

3860 

3 

3 

8.75 

5720 

7.  78 

75 

8.  87 

7200 

3a 

7.  77 
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8.89 
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3b 

5.  71 
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7.  29 

6220 

3c 

5.  70 

5680 

7.  28 
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4 

4 

9.  10 
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7.90 

7700 

9.34 
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5 

4 

7.  26 

4140 

5.  27 
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6.  95 

5590 

5a 

4 

7.  54 
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5.93 
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7. 04 

5660 

5b 

4 

7.00 
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5.  63 
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7.  26 

5850 

5c 

4 

6.85 

4350 

5.  59 
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7.  19 

5800 

5d 

4 

6.01 

4040 

5.  14 

4770 

6.84 

5510 

6 

- 

9.  50 

6110 

8.  32 

8170(3) 

9.70 

7920 

7a 

5 

4.  89 

3590 

4.  40 

4250 

5.  72 

4860 

7b 

5 

4.  53 

3340 

4.  27 

4140 

5.  58 

4770 

8a 

6 

5.  56 

5490(3) 

6.  88 

5870 

8b 

6 

6.  45 

6320{3) 

7.  59 

6360 

8c 

6 

6.  43 

4560 

6.  88 

4290 

8d 

6 

3.60 

2945 

4.  29 

3550 

Disregarding  structural  capability 

*2)  Ref.  Table  2 
(3) 

'  '  Engine  stops  "removed"  to  prevent  instability 
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Figure  3  Sissenwine  Wind  Profiles  For  1%  Probability  During  Winter 
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Figure  5  Avidyne  Wind  Profiles  For  1%  Probability  During  Winter 
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Figure  6  Marshall  SFC  Wind  Profiles  For  Saturn  Vehicle 
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Figure  7  Tail  Wind  Response  For  Various  Wind  Profiles  (Vehicle  A) 
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Figure  11  Regression  Curve  Of  Rigid  Body  Side  Wind  Response  (Vehicle  A) 
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Figure  15  Rigid  Body  Response  For  Design  Wind  Profiles  —  Vehicle 
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Figure  16  Rigid  Body  Response  For  Design  Wind  Profiles  —  Vehicle  C 
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Figure  17  Design  Wind  Response  ~  Vehicle  A 
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Figure  20  Design  Diagram  For  Wind  Shear  Response  —  Vehicle  A 
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Figure  21  Design  Diagram  For  Wind  Shear  Response  Vehicle  B 
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Figure  23  Shear  Versus  Shear  Length  For  Wind  Profiles  With 
Maximum  Wind  Velocity  Of  300  FPS 
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SHEAR  LENGTH  (1000  FT) 


Figure  24  Shear  Versus  Shear  Length  For  Wind  Profiles  With 
Maximum  Wind  Velocity  of  225  to  260  FPS 
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APPENDIX  A 

DEVELOPMENT  OF  DESIGN  WIND  PROFILE  SERIES 

The  design  diagrams  for  wind  shear  response  (Figs.  20,  21,  and  22)  logically 
should  be  obtained  for  the  mean  integrated  area  and  its  confidence  limits  at  the 
desired  level  of  confidence.  However,  as  long  as  the  required  statistics  are 
not  available,  a  method  of  estimation  must  be  provided.  This  is  accomplished 
by  means  of  four  design  profile  series.  Each  of  the  four  basic  profiles  has 
a  rational  basis,  permitting  it  to  be  expanded  into  a  series  with  maximum 
wind  velocity  as  the  argument.  Detailed  information  for  their  construction  is 
provided  in  this  Appendix. 

DESIGN  PROFILES  D1 

This  profile  series  has  two  parameters,  maximum  wind  velocity  and  shear. 
These  two  parameters  are  related  as  shown  in  Table  1  except  that  only  the 
1000-ft.  shear  is  used  and  expanded  to  3000-ft.  shear  length.  This  results  in 
a  large  incremental  wind  velocity  with  a  high  shear.  It  is  a  very  extreme 
condition  resulting  in  a  minimum  integrated  area.  The  parameters  are: 

Sp,  :  1000-ft  shear  (sec  )  associated  with  V  (Table  1) 

1/  tnzix 

V,  =  V  -3000S„  (fps),  Wind  velocity  at  (H  -3000)  ft. 

1  max  D  '  r  '  cr 

e  H  2 

A  =Jo  Cr  VdH  =  0.5V1Hcr+  1500  Vmax  (ft  /sec) 
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Design  Profile  Dj 

These  profiles  are  critical  for  minimum  area  or  Hcr  in  the  range  of  30,  000  ft. 
For  Hcr  =  30,  000  ft. 

A  *  1500  x  (10  V,  +  V  )  ft2/sec. 

1  max 

Typical  values  of  the  parameters  are  tabulated  in  Table  5  and  plotted  in 
Figure  26. 

DESIGN  PROFILES  D2 

Design  profiles  D2  have  only  one  parameter,  maximum  wind  velocity.  The 
basic  profile  simply  is  a  constant  shear  from  zero  ground  wind  to  maximum 
wind  velocity  at  critical  altitude.  These  profiles  are  critical  for  maximum 
integrated  area  and,  therefore,  the  critical  altitude  should  be  large,  40,  000  ft 

or  greater.  The  integrated  area  is  given  by: 

..  H  , 

A  a  f  cr  VdH  =  0. 5  V  H  (ft  /  sec) 

•  _  max  cr 


t 
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Design  Profile 

Typical  values  of  this  area  are  presented  in  Table  6. 


DESIGN  PROFILES  D3  AND  D4 


The  design  wind  profile  series  D3  and  D4  (see  Fig.  27  for  examples)  are 
synthetic  reconstructions  of  atmospheric  wind  profiles  using  the  AN/GMD-2 
data  from  Ref.  1  (Table  1).  They  correspond  to  the  Sissenwine  profile 
No.  3  (design  profile  D3,  300  is  identical)  and  differ  from  this  one  only 
in  critical  altitude  (D4  series)  and  the  property  to  have  a  specific  wind 
profile  rationally  related  to  any  maximum  wind  velocity.  The  purpose  of  the 
D3  and  D4  profile  series  is  to  approximate  the  mean  of  the  statistical  dis¬ 
tribution  of  the  integrated  area.  Both  profile  series  D3  and  D4  are  used  in 
order  to  obtain  a  "band"  which  may  be  expected  to  bracket  the  mean  integrated 
area.  Therefore,  the  D3  series  is  associated  with  smaller  integrated  areas 
(relatively  low  critical  altitude)  and  the  D4  series  with  larger  integrated 
areas  (relatively  high  critical  altitude).  Another  reason  for  the  two  series 
is  that  the  altitude  of  maximum  response  and  the  most  likely  altitude  of 
maximum  wind  velocity  usually  do  not  coincide.  With  two  series,  the  effect 
of  both  may  be  accounted  for. 
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o  V0  v,  v2  v3  V4  vMAX 

WIND  VELOCITY  (FPS) 


Design  Profile  and 

where  Vmax'  S1000*  S3000'  and  S5000  are  obtained  from  Table  l-  TYPical 
values  of  the  parameters  are  tat  netted  in  Table  7  and  plotted  in  Figure  28. 
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Table  5 

PARAMETERS  FOR  DESIGN  PROFILES  D1 


191  0.0280  107.0  1.892 


97.  7 
87.0 
73.  7 
58.  1 
41.5 


1.  720 
1.  530 
1.  300 
1.  037 
0.  755 
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Table  6 

PARAMETERS  FDR  DESIGN  PROFILES  D2 


V  fys 

mix" 

A  x  10‘6.  ft^/sec 

Hcr  =  40,000  ft 

H  =  45,  000  ft 

cr 

H  =  50,  000  ft 

300 

6,  00 

6.  750 

7.  500 

250 

5.00 

5.625 

6.250 

225 

4.  50 

5.063 

5.625 

191 

3.  82 

4.  302 

4.  775 

170 

3.40 

3.  825 

4.  250 

150 

3.  00 

3.  375 

3.  750 

130 

2.60 

2.925 

3.  250 

110 

2.  20 

2.475 

2.  750 

88 

1.  76 

1.980 

2.  200 
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PARAMETERS  FOR  DESIGN  PROFILES  D3  AND  D4 
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